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o ^ 

^ 1 Introduction 

o 

O ! In this talk, we report the results discussed in P, related to the phase structure of hot 

quark matter in presence of a background of chiral charge density, n^. The latter is 
introduced in the simplest way possible, namely by virtue of a chemical potential, /is, 
^ conjugated to ns. In more detail, after a brief introduction and list of motivations of 

p • this kind of study, we discuss the interplay between chiral symmetry restoration and 

p ^! deconfinement at finite /is, as well as the critical endpoint in the phase diagram and 

D \ its possible relationship with the critical endpoint of the phase diagram od Quantum 

Chromodynamics (QCD). In the talk, due to time limitation, we can emphasize few 
results related to the latter topic. Therefore we need to leave apart several applications 
^ ■ of the ideas, as well as of the formalism, developed here to the physics of heavy ion 

. collisions, with particular reference to the Chiral Magnetic Effect [2113111]. The latter 

^ \ might be relevant for the phenomenology of heavy ion collisions, because a copious 

\ production of gluon configurations, the QCD sphalerons, with a finite winding number 

is expected in the quark-gluon-plasma phase of QCD, see [5] and references therein. 
Because of the chiral Ward identity, the interaction of the sphalerons with the quarks 
causes a chirality change of the latters. As a consequence, a copious production of 
local domains in which chirality is imbalanced, is expected in the quark-gluon-plasma. 
\ The critical endpoint, CP, of QCD [6] is the cornerstone of the phase diagram of 

5^ ■ strongly interacting matter. At CP, a crossover line and a first order line are supposed 

to intercept. It is thus not surprising that an intense experimental activity is nowadays 
dedicated to the detection of such a point, which involves the large facilities at RHIC 
and LHC; moreover, further experiments are expected after the development of FAIR 
at GSI. Several theoretical signatures of CP have been suggested [HIH]. Despite the 
importance of CP, a firm theoretical evidence of its existence is still missing. In 
fact, the sign problem makes the Lattice Monte Carlo simulations difficult, if not 
impossible, in the large baryon-chemical potential (/x) region for iVc = 3 [9], see [10] 
for a recent review. Therefore, it has not yet been possible to prove unambiguously 
the existence and the location of CP starting from first principles simulations of 
grand-canonical ensembles. Moreover, the predictions of effective models are spread 
in the T — /i plane, see for example [HI 
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Interesting overcomings of the sign problem for the quest of CP are: analytic 
continuation of data obtained at imaginary chemical potential, fij [131 EH 113 116] : 
simulations at finite isospin chemical potential, see for example [171 UHl IlH]; sim- 
ulations in canonical, rather then grand-canonical, ensembles [20]; strong coupling 
expansion of Lattice QCD [211 122] • On the purely theoretical side, it has been sug- 
gested [23] that the use of orbifold equivalence in the large Nc approximation of QCD 
can lead to relations between the coordinates of CP at finite chemical potential, with 
those at finite isospin chemical potential. 

In this talk, we present the idea suggested in [1] about a new theoretical way to 
detect the CP, by means of Lattice simulations with Nc = 3. In order to accomplish 
this important program, we suggest to simulate QCD with a chiral chemical potential, 
/X5, conjugated to the chiral charge density, = tir — til, see [H [211 ISSl [261 121] for 
previous studies. Our idea, supported by concrete calculations within acmicroscopic 
effective model, is that CP can be continued to a critical endpoint at /is 7^ and /i = 0, 
that we denote by CP5, the latter being accessible to iV^ = 3 Lattice QCD simulations 
of grand-canonical ensembles [3l [28] . Therefore, the detection of the former endpoint 
via Lattice simulations, can be considered as a signal of the existence of the latter. 

The model used in the calculation, namely the Nambu-Jona-Lasinio model with 
the Polyakov loop [29] (PNJL model in the following) with tree level coupling among 
chiral condensate and Polyakov loop [30], gives numerical relations among the coor- 
dinates of CP5 and those of CP. In particular, the critical temperature turns out to 
be almost unaffected by the process of continuation; the critical value of the chemical 
potential, /Xc, on the other hand turns out to be almost half of the critical chiral 
chemical potential, n^c- 

Before discussing our results, it is important to spend some word more about the 
chiral chemical potential. In particular, we are aware that world at finite /is should 
be considered as a fictional one. As a matter of fact, /is cannot be considered as 
a true chemical potential because, in the confinement phase, the chiral condensate 
(gg) mixes left- and right-handed components of the quark field, leading to non- 
conservation of ns. Moreover, the quantum chiral anomaly leads to fiuctuations of 
the topological charge, which in turn causes the changes of the chiral density because 
of the Ward identity. Therefore, the point of view that we adopt is to consider /is 
as a mere mathematical artifice. However, the world at finite /is with A^^^ = 3 can 
be simulated on the Lattice. Therefore, it is worth to study it by grand-canonical 
ensemble simulations: it might furnish an evidence of the existence of the critical 
endpoint in the real world. 
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2 The Model with the Polyakov loop 



Because of its non-perturbative nature, we cannot make first principles calculations 
within QCD in the regimes to which we are interested in, namely moderate T, /z and 
/X5. Hence we need to rely on some effective model, which is built in order to respect 
(at least some of) the symmetries of the QCD action. To this end, we make use of 
the Nambu-Jona-Lasinio model [31] (see [32] for reviews) improved with the Polyakov 
loop [29], dubbed PNJL model, which has been used many times in recent years to 
describe successfully the thermodynamics of QCD with two and two-plus-one flavors, 
see [30l |33l [Ml [351 ISSl [371 [381 [Ml HO] and references therein. The model is interesting 
because it allows for a self-consistent description of spontaneous chiral symmetry 
breaking; even more, it allows for a simultaneous computation of quantities sensible 
to confinement and chiral symmetry breaking. We restrict here to a brief summary 
of the main equations; we refer to \T\ for a more detailed discussion. 

In the PNJL model, quark propagation in the medium is described by the following 
lagrangian density: 

C = q [i^D, - m + /X57°7' + m') Q + G [{qqf + {m^rq)'] ; (1) 

In the above equation, q corresponds to a quark field in the fundamental representa- 
tion of color group SU{3) and flavor group SU{2). We have a introduced chemical 
potential for the quark number density, /i, and a pseudo-chemical potential conju- 
gated to chirality imbalance, fi^. The chiral charge density, = ur — ul, represents 
the difference in densities of the right- and left-handed quarks. The imbalance of 
chiral density can be created by instanton/sphaleron transition in QCD, see [3] and 
references therein. At finite /is, a chirality imbalance is created, namely ns 7^ 0. For 
example, in the massless limit and at zero baryon chemical potential one has |3j 

n. = ^,+^. (2) 
37r^ 3 

If quark mass (bare or constituent) is taken into account, the relation n^{^^) cannot 
be found analytically in the general case, and a numerical investigation is needed, see 
for example p5]. 

In our computation we follow the idea implemented in [30j, which brings to a 
Polyakov-loop-dependent coupling constant: 

G = g[l-aiW -a2{L' + {L^f)\ , (3) 

The ansatz in the above equation was inspired by [121 |33] in which it was shown 
explicitly that the NJL vertex can be derived in the infrared limit of QCD, it has 
a non-local structure, and it acquires a non-trivial dependence on the phase of the 
Polyakov loop. We refer to [30] for a more detailed discussion. This idea has been 
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analyzed recently in [33], where the effect of the confinement order parameter on 
the four-fermion interactions and their renormalization-group fixed-point structure 
has been investigated. The numerical values of ai and 02 have been fixed in [30] 
by a best fit of the available Lattice data at zero and imaginary chemical potential 
of Ref. [^ 136]. In particular, the fitted data are the critical temperature at zero 
chemical potential, and the dependence of the Roberge- Weiss endpoint on the bare 
quark mass. The best fit procedure leads to ai = 0:2 = « = 0.2 ± 0.05. 

In the one-loop approximation, the effective potential of this model is given by 



V = U{L,L\T) + a^G- N^NfY^ 



=±1 



where 

Us = \J {\p\s - ii^f + ml , (5) 

corresponds to the pole of the quark propagator, and 

F_ = 1 + 3Le-'^('^=~^) + ?,L^e-^^^^^-^^ + e'^^^'^^-'') , (6) 
F+ = 1 + 3L^e-^("=+^) + 3Le-2'^("'^+'^) + e''^^^'^'^^'^ , (7) 

denote the statistical confining thermal contributions to the effective potential; Ug is 
given by Equation ([S]), with rriq = m — 2Ga. Once again the vacuum fluctuation term 
is regularized by means of a ultraviolet cutoff, that we denote by M. The relation 
between the chiral condensate and a in the PNJL model is a = (gg). 

We notice that the PNJL model considered here, which is dubbed Extended- 
PNJL in [20], has been tuned in order to reproduce quantitatively the Lattice QCD 
thermodynamics at zero and imaginary quark chemical potential. Hence, it represents 
a faithful description of QCD, in terms of collective degrees of freedom related to chiral 
symmetry breaking and deconfinement. 

The potential term U in Eq. (j3D is built by hand in order to reproduce the pure 
gluonic lattice data with Nc = ?> [33]. We adopt the following logarithmic form, 

W[L,I,T] = r^|-^lL + 6(r)ln[l-6lL + 4(l' + L=') -3(lL)2]| . (8) 
We refer to [33l [1] for the numerical values of the parameters used in this study. 



3 Critical endpoint at zero chemical potential 

In Figure [1] we plot the phase diagram of the model in the fi^ — T plane, for the 
case fi = 0. At any value of T and /X5, the chiral condensate and the Polyakov loop 
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Figure 1: [Color online). Phase diagram of the PNJL modeL The scale T^? = 173.9 
MeV corresponds to the critical temperature at /is = 0. 



expectation value are computed by a minimization procedure of the full potential dl]). 
The structure of our phase diagram is in agreement with previous model studies, 
see [2H |25]. Since chiral symmetry is broken explicitly by the quark mass and the 
phase transitions are actually crossovers, we identify the critical temperature with 
that at which dL / dT is maximum. We have checked that the latter deviates from that 
at which \da/dT\ is maximum only of a few MeV, in the whole range of parameters 
studied. With an abuse of nomenclature, we dub the pseudo-critical lines as second 
order and first order. It is clear from the context that the term second order transition 
has to be taken as a synonym of smooth crossover; similarly, the term first order 
transition is a synonym of discontinuous jump of the order parameters. 

In the Figure [H the grey dashed line corresponds to a smooth crossover. The solid 
line, on the other hand, denotes the first order transition. The dot corresponds to 
CPs- In the PNJL model we have access to the chiral condensate and to the Polyakov 
loop expectation value. As a consequence, we can label the phases of the model in 
terms both of confining properties, and of chiral symmetry. In the model at hand, 
because of the entanglement in Equation ([3]), the deconfinement and chiral symmetry 
restoration crossovers take place simultaneously. The region below the pseudo-critical 
line is characterized by confinement and spontaneous breaking of chiral symmetry; 
we label this phase as the confinement phase. On the other hand, the phase above 
the critical line is identified with the Quark-Gluon-Plasma phase. In this case, CPs is 
both chiral and deconfinement critical endpoint. For what concerns the coordinates 
of CPs we find, for the PNJL model. 
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Figure 2: {Color online). Evolution of the critical endpoint in the fj, — fi^ — T space, 
for the PNJL model. 



where T° = 173.9 MeV is the deconfinement temperature at /i = /^s = 0. 

Next we turn to discuss the more general case with both fi^ and fi different from 
zero. Our scope is to show that, at least within the models, CP naturally evolves 
into CP5. In particular the PNJL model, which is in quantitative agreement with the 
Lattice at zero chemical potential, gives a numerical relation among the coordinates 
of CP and CP5, which might be taken as a guide to estimate the coordinates of CP 
in QCD, once CP5 is detected. 

In Figure [2] we collect our data on the critical point of the phase diagram in the 
fi — fi^ — T space, in the case of the PNJL model. The orange solid line is the union 
of the critical points computed self-consistently at several values of fi: at any value 
of jj,, a point on the line corresponds to the critical point of the phase diagram in the 
/X5 — T plane. Thus the line pictorially describes the evolution of the critical point of 
the chiral model at hand, from CP to CP5. In the same Figure we plot a projection 
of the critical endpoint evolution curve onto the /i — /is plane, for the PNJL model. 
The indigo solid line corresponds to the /is-coordinate of the critical endpoint. The 
critical temperature is not so much affected when we continue CP5 to CP (we measure 
a change approximately equal to the 3%), therefore the projection in the /x — T plane 
is redundant. Our numerical results suggest the following relationships between the 
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critical coordinates at /i = and /is = 0: 

T 

— ^ 0.53 , ^ 0.97 , (PNJL) . (10) 

The model predictions (fTOj) relate the coordinates of CP to those of CP5. In 
particular, it is interesting that the critical temperature is almost unchanged in the 
continuation of CP to CP5. Of course, since these results are deduced by a model, it 
is extremely interesting and important to study how Equation (ITUI) is affected by the 
value of the bare quark mass, as well as by further interactions in the vector and axial- 
vector channels. These topics will be the subject of a forthcoming publication |47j. 
It is worth to anticipate some of the results of [17], namely that a larger value of the 
quark mass, as well as the vector interaction, move CP5 to larger values of fi^. The 
combination of these two factor, together with the finite size of the lattice cell, might 
explain the absence of CP5 in the Lattice simulations [28j. 



4 Conclusions and Outlook 

In this talk, we have reported on our results about the phase structure of hot quark 
matter in presence of a background of chiral charge density. Such a background is 
introduced by virtue of a chemical potential, yUs, conjugated to the chirality imbalance, 
'^5 = '^R—ni- Because of the fluctuations of the topological charge, which is connected 
with chirality imbalance in QCD via the quantum anomaly, fi^ should be treated as 
a pure mathematical artifice, and cannot be considered as a true chemical potential. 

This study is partly motivated by the potential applications to the Chiral Magnetic 
Effect [21 [3l H]. The latter might be relevant for the phenomenology of heavy ion 
collisions, because a copious production of gluon configurations, the QCD sphalerons, 
with a finite winding number is expected in the quark-gluon-plasma phase of QCD, 
see [S] and references therein. Because of the chiral Ward identity, the interaction 
of the sphalerons with the quarks causes a chirality change of the latters. As a 
consequence, a copious production of local domains in which chirality is imbalanced, 
is expected in the quark-gluon-plasma. 

After an overview on the phase diagram of hot quark matter at finite /is, obtained 
within an effective model, we have suggested the possibility of continuation of the 
critical endpoint of the phase diagram of A^^c = 3 QCD, CP, to a critical endpoint 
dubbed CP5 at finite /is and /i = 0. The worldsheet Ws = {/i = 0,/is = 0} has 
the merit that it can be simulated on the Lattice [3l [28] for Nc = 3. Even if Lattice 
results have been already published [28j, more care should be taken in the derivation 
of the relation between rzs and /is, since is a nonconserved quantity, hence it suffers 
renormalization effects which should be taken into account. 

The phase structure that we have discussed here is based on the PNJL model 
with entanglement vertex, introduced in [30], which offers a description of the QCD 
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thermodynamics in terms of collective degrees of freedom, which is in quantitative 
agreement with Lattice data at zero and imaginary chemical potential. 

One of our ideas is that simulations in the worldsheet W5 might reveal the ex- 
istence of a critical endpoint, CP5, in the phase diagram. Then, this critical point 
might be interpreted as the continuation of the critical point which is expected to 
belong to the phase diagram of real QCD, because of the continuity summarized in 
Fig. |2j Hence it would be an indirect evidence of the existence of the critical point in 
real QCD. 

In our calculations there are some factors that we have not included for simplicity, 
and that affect the location of CP5. For example, the bare quark mass and the vector 
interaction move CP 5 to higher values of n^. These observations might be helpful to 
understand why in the Lattice simulations of [28], no critical endpoint is detected. 
We plan to report on the aforementioned topics in the next future |17]. 
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